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ABSTRACT Certain short polycations, such as TAT and polyarginine, rapidly pass through the plasma 
membranes of mammalian cells by an unknown mechanism called transduction as well as by endocytosis 
and macropinocytosis. These cell-penetrating peptides (CPPs) promise to be medically useful when fused 
to biologically active peptides. I offer a simple model in which one or more CPPs and the phosphatidylser- 
ines of the inner leaflet form a kind of capacitor with a voltage in excess of 180 mV, high enough to create 
a molecular electropore. The model is consistent with an empirical upper limit on the cargo peptide of 
40-60 amino acids and with experimental data on how the transduction of a polyarginine-fluorophore into 
mouse C2C12 myoblasts depends on the number of arginines in the CPP and on the CPP concentration. 
The model makes three testable predictions. 



I. CELL-PENETRATING PEPTIDES 

In 1988, two groups [T] E] working on HIV reported 
that the trans-activating transcriptional activator (TAT) 
of HIV-1 can cross cell membranes. The engine driving 
this 86-aa cell-penetrating peptide (CPP) is its residues 
48-57 GRKKRRQRRR which carry a charge of +8e. Other 
CPPs soon were found. Antp (aka Penetratin, PEN) 
is residues 43-58 rqikiwfqnrrmkwkk of Antennape- 
dia, a homeodomain of the fly; it carries a charge of 
+7e. The polyarginine (Arg)„ carries charge +ne, where 
often n — 7, 8, or 9. Other CPPs have been discov- 
ered (VP22) or synthesized (transportan) . The struc- 
tural protein VP22 of the tegument of herpes simplex 
virus type 1 (HSV-1) has charge +15e. Transportan 
GWTLNS AG YLLG-K-INLKALAALAKKIL- amide is a chimeric 
peptide constructed from the 12 N-terminal residues of 
galanin in the N-terminus with the 14-residue sequence 
of mastoparan and a connecting lysine [3] . With its ter- 
minal amide group, its charge is +5e. 

These and other short, positively charged peptides can 
penetrate the plasma membranes of live cells and can tow 
along with them cargoes that greatly exceed the 600 Da 
restriction barrier. They are promising therapeutic tools 
when towing cleverly chosen peptide cargoes of from 8 to 
33 amino acids @] U El 13 HI 13 ESI EH E3 E3 13 ED3 EU ■ 

Many early experiments on CPPs were wrong be- 
cause the cells were fixed or insufficiently washed. Even 
careful experiments sometimes have yielded inconsistent 
results in part because fluorescence varies with the 
(sub)cellular conditions and the fluorophores [17] . 

Yet some clarity is emerging: TAT carries cargoes 
across cell membranes with high efficiency by at least two 
functionally distinct mechanisms according to whether 
the cargo is big or small [18] . Big cargoes, such as pro- 
teins or quantum dots, enter via caveolae endocytosis and 
macropinocytosis [THJ ED] > and relatively few escape the 
cytoplasmic vesicles in which they then are trapped |18j . 

Small cargoes, such as peptides of fewer than 30- 
40 amino acids, enter both slowly by endocytosis and 
rapidly by transduction with direct access to the cytosol, 



an unknown mechanism that uses the membrane poten- 
tial pi] El E2 EH El]. Peptides fused to TAT enter cells 
within seconds [25 . 

It remains unclear how big cargoes aided by several 
CPPs enter cells EE]. For instance, superparamegnetic 
nanoparticles encased in aminated dextran and attached 
to 45 tat peptides are thought to enter cells by adsorptivc 
endocytosis (27] EH] ES] but they do enter slowly at 4° 

c [U]. 

This paper is exclusively about how polycationic cell- 
penetrating peptides, specifically oligoarginines, trans- 
duce small cargoes directly into the cytosol. Sec.[TT]recalls 



III 



some basic facts about plasma membranes, and Sec 
explains why ions do not normally pass through plasma 
membranes. Sec. |IV| describes a simple model of the 
transduction of CPPs in which electroporation and phos- 
phatidylserine play key roles. In this model, one or more 
positively charged CPPs on the outer leaflet and the neg- 
atively charged PSs under it on the inner leaflet form a 
kind of capacitor, which enhances the membrane poten- 
tial to a voltage in excess of 180 mV, which is sufficient 
to create an electropore. Sec. [V] shows that the model 
is consistent with an empirical upper limit on the cargo 
of 40-60 amino acids and with measurements made by 
Tiinncmann et al. |31j on the fraction of mouse myoblasts 
transduced by polyarginines carrying fluorophores of 400 
Da. Sec. |VI| tells how to test three predictions of the 
model. The paper ends with a short summary in Sec. |VIl[ 



II. MAMMALIAN PLASMA MEMBRANES 

The plasma membrane of a mammalian cell is a lipid 
bilayer that is 4 or 5 nm thick. Of the four main phos- 
pholipids in it, three — phosphatidylethanolamine (PE), 
phosphatidylcholine (PC), and sphingomyelin (SM) — are 
neutral, and one, phosphatidylserine (PS), is negatively 
charged. In live cells, PE and PS are mostly in the cy- 
tosolic layer, and PC and SM in the outer layer [32] [33] . 
Aminophospholipid translocase (flippase) moves PE and 
PS to the inner layer; floppase slowly moves all phospho- 
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lipids to the outer layer [32"] . 

Incidentally, the surfaces of bacteria are different. The 
cell wall of a Gram-positive bacterium (e.g., Streptococ- 
cus or Staphylococcus) is covered with negatively charged 
teichoic acids; the outer leaflet of the outer membrane of 
a Gram-negative bacterium (e.g., E. coli or Salmonel- 
lum) is tiled by negatively charged lipopolysaccharidcs 
(LPS) held together by divalent cations [34"1 135] . 

Glycolipids make up about 5% of the lipid molecules 
of the outer layer of a mammalian plasma membrane 
where they may form lipid rafts. Their hydrocarbon 
tails normally are saturated. Instead of a modified phos- 
phate group, they are decorated with galactose, glucose, 
GalNAc = N-acetylgalactosamine, and other sugars. The 
most complex glycolipids — the gangliosides — have nega- 
tively charged sialic-acid (NANA) groups. Incidentally, 
cholera toxin binds to and enters cells that display the 
Gmi ganglioside. [33^ 

A living cell maintains an electrostatic potential of be- 
tween 20 and 120 mV across its plasma membrane. The 
electric field E within the membrane points into the cell 
and is huge, about 15 mV/nm or 1.5 x 10 7 V/m if the 
potential difference is 60 mV across a membrane of 4 
nm. Conventionally, one reports membrane potentials 
as the electric potential inside the cell minus that out- 
side, so that here AV = —60 mV. Near but outside 
the membrane, this electric field falls-off exponentially 
E(r) = E cxp(—r/D() with the ratio of the distance r 
from the membrane to the Debye length Dg, which is of 
the order of a nanometer. The rapid entry of TAT fused 
to peptides is frustrated only by agents that destroy the 
electric field E [TB], which applies a force qE to a CPP 
of charge q. 

Most of the phospholipids of the outer leaflet of the 
plasma membrane are neutral PCs & SMs. They vastly 
outnumber the negatively charged gangliosides, which are 
a subset of the glycolipids, which themselves amount only 
to 5% of the outer layer. Imagine now that CPP-cargo 
molecules are in the extra-cellular environment. Many 
of them will be pinned down by the electric field E(r) 
just outside the membrane, their positively charged side- 
chains interacting with the negative phosphate groups of 
neutral dipolar PC & SM head groups [23]. (Other CPP- 
cargo molecules will stick to negatively charged ganglio- 
sides and to glycosaminoglycans (GAGs) attached to 
transmembrane proteoglycans (PGs); these slowly will 
be endocytosed. PGs with heparan-sulfate GAGs are 
needed for TAT-protein endocytosis [35] ■ A more de- 
tailed analysis than that of this work might model the 
effect of these anionic matrix compounds upon transduc- 
tion.) It is crucial that the dipolar PC & SM head groups 
are neutral and so do not cancel or reduce the positive 
electric charge of a CPP-cargo molecule. The net posi- 
tive charge of a CPP-cargo molecule and the negatively 
charged PSs under it on the inner leaflet form a kind of 
capacitor. This is the starting point for the model de- 
scribed in Sec. II VI 



III. THE PROBLEM 

The dielectric constant q ~ 2 of the hydrocarbons of 
a lipid bilayer is much less than that of water e w ~ 80. 
Thus, the difference AE w ^e in the electrostatic energy 
of an ion of charge q and effective radius a in the bilayer 
and in water [37] is 



AK, 



1 



1 



(1) 



or 3.5 eV if the ion's charge is that of the proton and 
its radius is a = 1 A. This energy barrier is far larger 
than the 0.06 eV gained when a unit charge crosses a 60 
mV phospholipid bilayer. Thus, an ion will not cross a 
cell's plasma membrane unless a transporter or a channel 
facilitates (and regulates) its passage. 

In the present model of CPP transduction, the electro- 
statics of the CPP-cargo complex and the role of PSs on 
the inner leaflet play key roles. 

The electrostatics of a cationic polypeptide such as 
TAT or polyarginine are more complex than for an ion. 
I will model the CPP and its cargo in water as a sphere 
with its positive charges on its surface. The density of a 
protein of mass M kDa is estimated [38] to be 



p(M) = (fj.8491 + 0.0873 e~ M/13 ) kDa/r 



(2) 



A CPP-cargo complex would not be expected to fold as 
densely as a natural globular protein, and so for it the 
estimate p(M) is something of an upper bound. The 
radius r of a putative sphere consisting of M kDa of 
CPP and cargo then would be 



> 



3 M \ 1/3 
4^ p(M) ) 



(3) 



For instance, a CPP of N arginines and a tiny fluorophorc 
cargo of 400 Da has a mass of M N = 0.1562iV + 0.4 kDa, 
and so its radius would satisfy 



> 



M N 



An p(M N ) 



1/3 



(4) 



or r = 0.75 nm for N — 8 arginines. The lower bounds 
on the radii for N = 5-12 are listed in column 2 of the 
Table ffl). 

For larger cargoes of A = 50-100 amino acids of 130 
Da each, the lower bounds on the radii range from 1.25 to 
1.59 nm. (In what follows, A will represent the number 
of amino acids in the cargo or the mass of the cargo in 
Daltons divided by 130 Da.) Adding another 0.8 nm 
for the PC/SM head groups would extend these lower 
bounds on the radii to 2.05-2.39 nm. 

If the CPP-cargo molecule were a charged conducting 
sphere of radius r and charge q, then its electrostatic 
energy in water would be 



E(N,A,q,w) 



(5) 
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This term neglects the short-distance detail of the electric 
field near the q/e positive unit charges e of the CPP-cargo 
molecule. So a short-distance correction term 



E sdc (a,q,w) 



qe 



8iT€Qe w ci 



(6) 



proportional to q must be added to E(N,A,q,w). The 
short distance a is a parameter, which will turn out to be 
a few A because the term E s d c is a correction to be added 
to E(N, A, q, w) and not the entire electrostatic energy. 

The electrostatic field of the cell attracts the CPP- 
cargo molecule to the surface of the cell. While on the 
outer leaflet of the plasma membrane, the electrostatic 
energy of the CPP-cargo molecule and its short-distance 
correction are no longer given by their values in water 
Eqs. ([5] & [6| but instead are those appropriate to the 
interface between water and lipid 



E(N, A, q, wi) = 



87reoer 



and 



E sdc (a,q,w£) 



qe 



87re ea 



where e is the mean permittivity 
1 < 

e = 7i I 6 ™ + e e) 



(7) 



(8) 



(9) 



The CPP-cargo molecule enters the lipid bilayer as a 
CPP-cargo-PC/SM complex with the phosphate groups 
of the PC and SM of the outer leaflet bound to the 
positively charged guanidinium and amine groups of the 
CPP [23]. The positive charges of the phosphocholinc 
groups of PC and SM are about d = 5 A from their 
phosphate groups [39] ■ The binding of PC and SM there- 
fore approximately increases the effective radius of the 
charged sphere to r m as r + d. The electrostatic energy 
of this complex in the hydrocarbon tails of the lipid bi- 
layer then is 



E(N,A,q,£) 



87re Q(Y + d) 
apart from a short-distance correction factor 



E sdc (a,q,£) = 



qe 



87reoQa 



(10) 



(11) 



similar to ([6]). 

Apart from correction terms, the electrostatic energy 
penalty when the CPP-cargo molecule enters the lipid 
bilayer from water as a CPP-cargo-PC/SM complex is 
the difference 

AE° w AN 7 A,q) « E(N,A,q,£)-E(N,A,q,w£) (12) 

r + dei" 



lire ec(r + d) 



1 - 



r e 



Because the thickness £ = 4.4 nm of the lipid bilayer 
is at most a few times the diameter of the CPP-cargo- 
PC/SM complex, we also must include the Parsegian cor- 
rection [37] 



AE P = q —ln' 2e * 



(13) 



which holds when a uniformly charged sphere is inserted 
into the middle of a lipid layer of thickness I. The sum of 
the water-to-lipid energy ( 12 1 and Parsegian's correction 
(fl3l is 



AE w/ (N, A, q) = AE° W /N, A, q) + AE f 



(14) 



The energy AE W ^(N, A, q) is listed in column 3 of Ta- 
ble for a CPP of N — 5-12 arginines towing a fluo- 
rophore cargo of 400 Da with d = 0.5 nm. 

The short-distance correction terms augment this 
penalty by 



AE sdc (a,q) = E sdc (a,q,£) - E sdc (a,q,w) 
8Tre e e a \ e w 



(15) 



and do not require Parsegian's correction because they 
are short-distance effects. This short-distance correction 
AE sdc is listed in column 4 of Table ^ for CPPs of N = 
5-12 arginines and a representative value of a = 4.5 A 
for the short-distance parameter. 

The net electrostatic energy penalty when the CPP- 
cargo molecule enters the lipid bilayer from water as a 
CPP-cargo-PC/SM complex is then the sum of ([12] [l3~ 



AE„ 



AE„, 



AE P + AEi 



dc ■ 



(16) 



A CPP of 8 arginines carrying a fluorophore of 400 Da 
(A = 3 w 400/130) has a radius r of 0.75 nm, and with 
a = 4.5 A, the change ( 16 1 in its electrostatic energy on 



going from water to lipid is 



AE„ 



16.9 eV. 



(17) 



This energy barrier is 35 times bigger than the energy 
0.48 eV that it gains by crossing a potential difference of 
60 mV. So how and why does it cross? 



IV. THE MODEL 

My answer is that one (or more) oligoarginines and 
the phosphatidylserines (PSs) of the inner leaflet together 
with their counterions form a kind of capacitor with an 
electric field strong enough to form a reversible pore in 
the plasma membrane. The transmembrane potential is 
the sum of three terms — the resting transmembrane po- 
tential AV ce u of the cell in the absence of CPPs, the 
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TABLE I: The radius r of a CPP-cargo molecule of N 
arginines and a cargo of 400 Da, its change in electrostatic 
energy AE W ^ when transferred from water to hydrocarbon 



(Eq. 141, and the short-distance correction AEsdc (Eq. 15 1 



Distances are in nm and energies in eV. 



N 


r 


AE W) i 


AE sdc 


5 


0.67 


4.61 


3.90 


6 


0.70 


6.39 


4.68 


7 


0.73 


8.41 


5.46 


8 


0.75 


10.64 


6.24 


9 


0.78 


13.06 


7.02 


10 


0.80 


15.68 


7.80 


11 


0.82 


18.48 


8.58 


12 


0.84 


21.44 


9.36 



transmembrane potential AVcpp due to an oligoargi- 
nine, and the transmembrane potential AVnuCI due to 
the countcrions of the extracellular medium 



AV = AV, 



cell 



AVr 



CPP 



AVr 



NaCh 



(18) 



The resting transmembrane potential AV ce u of the cell 
varies between about 20 mV to more than 70 mV, de- 
pending upon the type of cell. Ideally, it is measured 
experimentally. 

My model is based upon several considerations, which 
I discuss in turn in this section. The first subsection 
describes the basic facts about electroporation. The sec- 
ond subsection presents the electric potential V due to 
a charge in the extracellular medium; the derivation of 
that potential is in an appendix. This potential im- 
plies that charges on opposite sides of the lipid bilayer 
are effectively decoupled, which simplifies the subsequent 
analysis. The third subsection describes a Monte Carlo 
simulation of the response of the phosphatidylserines of 
the inner leaflet to an oligoarginine interacting with the 
phosphate groups of the outer leaflet. To a very good 
approximation, the PSs are distributed uniformly and 
randomly because they are nearly decoupled from the 
oligoarginine. The fourth subsection uses the potential 
V to compute the contribution AVcpp of an oligoargi- 
nine to the transmembrane potential. The fifth subsec- 
tion describes a Monte Carlo simulation of the effect 
AVjvac; of the sodium and chloride ions in the extra- 
cellular medium upon the transmembrane potential, ex- 
tracellular medium. The section ends with a summary of 
the model. 



A. Electroporation 

Electroporation is the formation of pores in membranes 
by an electric field. Depending on the duration of the 
field and the type of cell, an electric potential difference 
across a cell's plasma membrane in excess of 150 to 200 
mV will create pores. There are two main components to 
the energy of a pore. The first is the line energy 27rr7 due 



to the linear tension 7, which is of the order of 10 -11 J/m. 
The second is the electrical energy — 0.5AC7rr 2 (AV) 2 in 
which Ay is the voltage across the membrane and AC = 
C w — Ct is the difference between the specific capacity 
per unit area C w — e w eo/t of the water-filled pore and 
that Ci = ef,€o/t of the pore-free membrane of thickness 
t. There also is a small term due to the surface tension 
£ of the plasma membrane of the cell, but this term 
usually is negligible since S is of the order of 2.5 x 10~ 6 
J /m 2 |40j . The energy of the pore in a plasma membrane 
is then gH HI H3 H US 



E(r) — 27TT7 - 



7rr 2 E 



\-nr 2 AC (AV) 2 



This energy has a maximum of 



E(r c ) 



7T7 



2tt7 2 



Y,+ \AC(AV) 2 AC(AV) 2 



at the critical radius 



AC (AV) 2 



27 

AC (A^) 2 ' 



(19) 



(20) 



(21) 



In Fig. [IJ the Boltzmann factor e -E{r)/(kT) ( xl0 o) is 
plotted as a function of the radius r of the pore up to 
r c for various transmembrane voltages from —200 (solid, 
red) to —400 mV (dot-dash, cyan). Clearly, the chance 
of a pore forming rises steeply with the magnitude of the 
voltage and falls with the radius of the pore. 

If the transmembrane potential AV is turned off before 
the radius of the pore reaches r c , then the radius r of the 
pore usually shrinks quickly (well within 1 ms [42]) to a 
radius so small as to virtually shut-down the conductivity 
of the pore. This rapid closure occurs because in ( 19 ) the 
energy 2nrj dominates over — 7rr 2 E, the surface tension 
£ being negligible. Such a pore is said to be reversible. 
But if AV remains on when r exceeds the critical radius 
r c , then the pore usually will grow and lyse the cell; such 
a pore is said to be irreversible. 



The formula (21 ) provides an upper limit on the radius 
of a reversible pore. This upper limit drops with the 
square of the transmembrane voltage Ay from r c = 3.6 
nm for AV = -200 mV, to 1.6 nm for Ay = -300, and 
to 0.9 nm for Ay = -400 mV. 

The time tg for a pore's radius to reach the critical 
radius r c is the time to lysis; it varies greatly and ap- 
parently randomly even within cells of a given kind. In 
erythrocytes, its mean value drops by nearly an order of 
magnitude with each increase of 100 mV in the trans- 
membrane potential |42j and is about a fifth of a second 
when Ay = -300 mV. 

In the present model, however, the potential is imposed 
by the CPP and the PSs, and so when that potential 
causes a pore to form, the CPP and its cargo may enter 
the cell through the pore that they have formed, arid once 
they do, the potential drops to its normal resting value, 
usually less than -100 mV, and the pore virtually closes 
within 1 ms. 
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1: The Boltzmann factor e 



-E(r)/kT 



(xlOO) for energy 
(Eq.(19l) is plotted against the radius r of the pore 



FIG. 

E(r) 

up to the critical radius r c (Eq.(21l) for the transmembrane 
voltage AV = -200 (solid red), -250 (dashes green), -300 
(dots blue), —350 (dots magenta), and —400 mV (dash-dot 
cyan). 



The oligoarginine(s) on the interface between the outer 
leaflet and the extra-cellular environment, the negatively 
charged head groups of the PSs below them in inner 
leaflet, and their counterions create an electric field and 
a transmembrane potential AV. The chance of this po- 
tential forming a pore of radius r is proportional to the 
Boltzmann factor e~ E ^^ kT \ which is plotted in Fig. [l] 
The higher the potential AV and the narrower the pore, 
the greater the chance of pore formation. 



B. The Potential of an External Charge 

As shown in Appendix A, the electrostatic potential in 
the lipid bilayer Vg(p,z) due to a charge q at the point 
(0, 0, h) on the z-axis a height h above the interface be- 
tween the lipid bilayer and the extra-cellular environment 



Vt(p,z) 



q 



47rene 



n=0 



a/,9 2 + {z + 2{n + fjt + h) 



^Jp 2 + {z~ 2nt - h) 2 
(22) 



for < h < t, in which t is the thickness of the lipid bi- 
layer, e w £ = (e w + eg) /2 is average of relative permittivity 
of the extra-cellular fluid e w and that of the lipid bilayer 



FIG. 2: The electric potential V(p, z 
for p = 1 nm as a function of the height z 



from |22j|25|l in Volts 
(nnij 



above the 

phospholipid bilayer (and in the extracellular environment) 
for a unit charge q = ej at (p,z) = (0,0) (top curve), (0,1) 
(middle curve), and (0,2) nm (bottom curve). The lipid bi- 
layer extends from z — to z = —5 nm, and the cytosol lies 
below 2 = — 5 nm. The relative permittivities were taken to 
be e w = e c = 80 and ei — 2. 



0.07 
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o 

Oh 
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2 4 6 8 10 
Distance p (nm) from z-axis 



FIG. 3: The electric potential V(p,z) from (24) in Volts for 
< p < 10 nm at a height z = 0.5 (nm) above the phospho- 
lipid bilayer for a unit charge q = ej at (p, z) = (0,0) nm. 
Same geometry and parameters as in Fig. [2] 
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q, and p and p' are the ratios 



P - 



and p' = 



eg 



(23) 



which lie between and 1. The potential in the extra- 
cellular medium is 



V w (p,z) 



1 



P 



4ire e w \r + (z + h) 2 



OO 

E 



p n-X p ln 



t ^ y/p 2 + (z + 2nt + hf 



(24) 



in which r = \J p 2 + (z — h) 2 is the distance from the 
charge q. The potential in the cytosol due to the same 
charge q is 



V c (p,z) 



qti 



-E 



(pp'Y 



4ne e we e ec ^ Q y^5 + ( z - 2nt - h) 2 



where tg c is the mean relative permittivity 
The first 1000 terms of the series (22 1, p4| 



(25) 



& (251 for 



the potentials Vi(p, z), V w {p,z), and V c (p, z) are plot- 
ted in Fig. [2] (in Volts) for p = 1 nm as a function of 
the height z (nm) above the phospholipid bilayer (and in 
the extracellular environment) for a unit charge q = |e| 
at (p,z) = (0,0) (top curve), (0,1) (middle curve), and 
(0,2) nm (bottom curve). The lipid bilayer extends 
from z = to z = —5 nm, and the cytosol lies below 
z = —5 nm. The relative permittivities were taken to 
be e w = € c = 80 and eg = 2. Fig. [3] plots the potential 
V w (p,z) in the extracellular region due to a unit charge 
at the origin as a function of p for z = 0.5 nm. 

In and near the extracellular region, these potentials 
are fairly well approximated by the simple formulas 



V w (p,z) 
Vi(p,z) 



1 



P 



Aire e w \ r ^p 2 + (z + h) 2 
g 

4ire e w i r 



(26) 
(27) 



which hold when the lipid bilayer is infinitely thick. But 



the potential drops significantly below this formula (27) 
as z descends deeper into the bilayer until it nearly van- 
ishes at the lipid-cytosol interface and in the cytosol. In 
fact, a charge of 12|e| at the origin raises the potential 
V^(p, z) on the interface at (p, z) = (1,-5) nm only to 
0.0079 V. Thus the energy advantage of a PS at (1, -5) 
nm is only 0.0079 eV, which is much less than kT^ ss 0.027 
eV. So a CPP on the interface between the lipid bilayer 
and the extracellular fluid has a very small effect on the 
PSs of the inner leaflet whose negative charges lie on the 
lipid-cytosol interface. 

It follows that the counterfoils of the extracellular fluid 
also have little effect upon the PSs. And since the electric 
permittivities of the extracellular fluid and of the cytosol 




-25 -20 -15 -10 -5 5 10 15 20 25 



FIG. 4: Superposition of 10 plots of the positions (x, y, —5) 
nm of 255 phosphatidylserines inside a disk of radius 25 nm 
on the lipid-cytosol interface 5 nm below an a-helix of 12 
arginines at (0, y, 0) for —0.96 < y < 0.8 nm (x's, green). To 
a good approximation, the PSs are uniformly and randomly 
distributed. 



are similar, we may view Fig.[2]upside-down and conclude 
that the PSs and the K + and Cl~ ions of the cytosol 
have little effect upon the CPP and the counterions of 
the extracellular fluid except to contribute most of the 
transmembrane potential that exists in the absence of 
CPPs. Charges in the cytosol are effectively decoupled 
from those in the extracellular environment. 

We may draw a further lesson from the sharp drop in Vi 
across the lipid bilayer shown in Fig. [2] The transmem- 
brane potential due to a CPP on the interface (z = 0) is 
much larger than the simple formula ( |27| would imply. 
This is why CPPs are transduced. 



C. Monte Carlo of the Phosphatidylserines 

Phosphatidylserines (PSs) make up some 8-18% of 
the inner leaflet by weight |46J . They diffuse laterally 
within that leaflet with a diffusion constant D 1=3 10~ 8 
cm 2 /sec [47] and so within one second spread to an area 
of 12 pm 2 , which is a significant fraction of the surface 
area of a eukaryotic cell. 

My Monte Carlo simulations of the distribution of the 
PSs of the inner leaflet verified the conclusions of the last 
subsection (IV B I based upon the analytic potentials (|22T 



25 1 and showed that the PSs are randomly and uniformly 
distributed, at least to a good approximation. Figure [4] 
superposes 10 snapshots of the (x, y) locations of 255 PSs 
in a disk of radius 25 nm that is 5 nm directly below a 
12-mer of arginine R 12 . The snapshots were taken every 
2000 sweeps after 25,000 thermalizing sweeps. The 255 
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PSs moved so as to minimize their free energy due to 
interactions with the R 12 , with each other, and with the 
PSs outside the disk. The distribution shows no obvious 
clustering. 

I considered the case of a single CPP of 5 < N < 
12 arginines. A CPP of N arginines can form an a- 
helix of length L a ss 0.16 (N — 1) nm, a random coil of 
length L r w 0.25 (N — 1) nm, or a /3-strand of length 
Lp m 0.34 (N — 1) nm. The random-coil and /?-strand 
configurations spread the positively charged guanidinium 
groups farther apart and so would be expected to cluster 
the PSs even less than the a-helix configuration. So in 
the simulations of this subsection, I only used the a-helix 
configuration. 

The Monte Carlo code [15] assumes that a PS has a 
cross-sectional area of 1 nm 2 and that the PSs make up 
13% of phospholipids of the inner leaflet. There are then 
about 0.13 ttR 2 PSs in a disk of radius R nm, or 255 PSs 
in a disk of radius 25 nm. The codes allow these 255 PSs 
to move about within that disk attracted by the electric 
potential of the N or 2N positively charged arginines 
and repelled by each other and by the PSs outside the 
disk, which are treated as a uniform surface charge. The 
computations are facilitated somewhat by the continuity 



of the electric potentials (22 24 & 25 1 across the inter- 
faces at z — and z — — t between the lipid bilayer and 
respectively the extra-cellular medium and the cytosol. 

The code assumes that the N arginines form an alpha 
helix with positive charges at the points 



r jC pp = (0, 0.16 (N/2 -j), 0). 



(28) 



The PSs were allowed to move in two dimensions within 
the disk of radius 25 nm in the inner leaflet at sites r k 
for k= 1,...2N. 

The electrostatic energy of a single PS at the point r k 
is the sum of three different energies 



E k = E, 



k,CPP 



E, 



k,PSs 



E, 



fe,<7* 



(29) 



The first energy Ek^cpp is that due to its interaction 
with arginines of the CPP(s) 



E, 



k.CPP 



M 

£ 

3=1 



eV c (\r k -r jC pp\,-t) (30) 



in which M = N for a single CPP and M = 2N for 
two CPPs. The second energy Ek.pss is that due to the 
interaction of the fcth PS with the Nps — 1 = 254 other 
PSs in the disk 



N PS 

E k ,pss= ^2 eV w (\r k -r k /\,Q). 

k'—l,k'^k 



(31) 



The third energy E ki(T is that due to the interaction of the 
fcth PS with all the PSs outside the disk represented by a 
uniform surface charge a = 0.13 e nm -2 . In Appendix B, 
I derive the approximation 



Eh, 



aR 
2e e v 



oo 
n—1 



l 



(2n) 



( n l)2 2 2r 



2 1 1 



(32) 



apart from an irrelevant infinite constant. In the com- 
puter programs, the upper limit on the summation was 
n = 800. 

The Monte Carlo codes use a simple Metropolis step 
in which the x-y coordinates of a single PS, the /cth, are 
randomly varied by as much as ±5 nm (to keep the ac- 
ceptance rate down to 68%). The codes accept any move 
that lowers the energy E k as given by (29 1 and also ac- 
cept any move that raises E k by AE k conditionally with 
probability 



p = e -AE k /(k B T) 



(33) 



in which kg is Boltzmann's constant and T is 37 Celsius. 
A sweep consists of k — 1 , . . . , Nps = 255 Metropolis 
steps. Each simulation started from a random configura- 
tion of PSs in the disk of radius R = 25 nm and ran for 
45,000 sweeps. Measurements began after 25,000 sweeps 
for thermalization. 

The radius of an electropore is about r p = 1 nm, and 
the thickness of the plasma membrane was taken to be 
t = 5 nm. The code measured the electrostatic potential 
across the lipid bilayer between points that were offset 
in the x-direction by 1 nm, that is, between the points 
(1,0,0) and (1,0,-5) nm. The code measured the volt- 
age across the membrane every 10 sweeps and recorded 
the positions of the PSs every 2,000 sweeps. 

The transmembrane potential V due to a single a-helix 
oligoarginine TL N and its cloud of PSs rises with the num- 
ber of arginines from about V — —380 mV for N = 5 
arginines to V — -630 mV for N = 9 and V = -800 
mV for N = 12. These voltages are so high that one 
would have expected electroporation even for R 5 s, which 
is not seen. The PSs contributed only about 70 mV to 
the transmembrane potentials listed in Table |VI| Thus, 
although they facilitate transduction, they are not re- 
sponsible for the very high voltages listed in the table. 
These voltages are too high because the simulations did 
not include the Na + and CI™ ions in the extra-cellular 
medium. 

In all these simulations, the mean value of the dis- 
tance of the PSs from the point (0, 0, -5) nm was about 
17 nm. The distributions of the PSs across the disk of 
radius r = 25 nm appeared uniform and random, with 
little clustering under the CPPs as shown in Fig. [4] The 
reason for the tepid PS response to the electric field of 
the CPPs can be seen in Fig [2] the electric potential 
V(p, z) drops off sharply as z descends through the lipid 
bilayer and is very small near the lipid-cytosol interface. 
This uniformity of the PS distribution on the inner leaflet 
means that we need not simulate their behavior explic- 
itly. We can use the resting transmembrane potential in 
the absence of CPPs to represent both the PSs and the 
counterions of the cytosol. Ideally, one should take it 
from experimental measurements. 
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TABLE II: The voltage differences AVcpp (mV) across the 
plasma membrane due to an oligoarginine as an a-helix, a 
random coil, or a /3-strand. Neither the salt potential AVno.ci 
nor the resting cell potential AV ce ii is included. 



N 


a-helix 


R^ random coil 


/3-strand 


5 


- 312 


- 302 


- 291 


6 


- 376 


- 362 


- 346 


7 


-439 


- 419 


- 393 


8 


- 502 


- 472 


- 425 


9 


- 562 


- 521 


-455 


10 


- 620 


- 557 


-476 


11 


- 676 


- 587 


-499 


12 


- 729 


- 614 


- 516 



D. The Potential of an Oligoarginine 

This subsection computes the transmembrane poten- 
tial AVcpp due to an R™ oligoarginine whose n unit 
positive charges for 5 < n < 12 were fixed at the points 



Tjcpp = (0, (N/2 - j) Ay, 0) 



(34) 



in which Ay = 0.16, 0.25, and 0.34 nm respectively for 
an a-helix, a random coil, and a /3-strand. I took this 
AVcpp to be the difference 



AV ( 



CPP 



= (V c (p,-t))-(V w (p,0)) 



(35) 



in which (V w (p, 0)) and (V c (p, —t)) are the mean values of 
the R"'s electric potential on two disks of radius r p = 1 
nm at z = and at z = —t. 

I used a Monte Carlo code [48 j to numerically integrate 
the appropriate potential V w or V c (Eqs. (24 or|25|) over 
the two disks. The code used a million random points on 
each of the disks (of which the fraction (4 — 7r)/4 = 0.215 
were discarded because they lay outside the d isk) . In 
this code, I kept 100 terms in the series p4| fc [25|; the 
error introduced by this truncation is completely negligi- 
ble (about 2 parts in 10 million). 

The resulting transmembrane potentials AVcpp are 
listed in Table [IT] The magnitude of AVcpp naturally 
increases with the charge n|e|. Because the n charges are 
more spread out in a /3-strand than in a random coil, the 
magnitude of AVcpp is less for a /3-strand than for a 
random coil of the same charge, and similarly for a coil 
and an a-hclix. 



E. Monte Carlo of the Counterions 

In this subsection, I use Monte Carlo methods to com- 
pute the transmembrane potential AVjvaCZ due to the 
sodium and chloride ions of the extracellular medium 
near an oligoarginine. 

The Na + , K+, Mg++, Ca ++ , and Cl~ concentrations 
in the extracellular medium respectively are 145, 5, 1-2, 
1-2, and 110 mM [15]. I approximated their effects by 



setting the Na + and Cl~ concentrations to 156 mM and 
ignoring the other ions. I used an active volume that was 
10 nm wide and 20 nm long, and that rose from the lipid 
bilayer to a height of 5 nm. In this active volume of 200 
(nm) 3 , I put 94 sodium ions and (94+n) chloride ions so 
as to make the charge within the active volume neutral. 

To prevent the sodium and chloride ions from avoiding 
the walls and ceiling of the active volume, I surrounded 
the walls and ceiling of the active volume with a 1000 
(nm) 3 5 nm-thick passive volume in which I randomly 
placed 470 Na+ and 470 Cl~ ions. 

The Monte Carlo code [15] used the potential V w of 



Eq. ( 24 1 to compute the energy of an individual sodium 
or chloride ion in the active volume due to its interaction 
with all the ions in the active and passive volumes and 
with the CPP(s) which did not move. The fixed positions 
(34 1 of the n charges of the oligoarginine depended upon 



whether the R n was configured as an a-helix, a random 
coil, or a /3-strand. The ions in the passive volume also 
didn't move, retaining their original random positions, 
which were different in each run. To speed up the com- 
putation, I used only the first 8 terms in the series (24 1 
for V w (p, z) , which introduced an error of about 0.6%"^ 
In order to prevent the Na + and Cl _ ions from collaps- 
ing into neutral composite particles of infinite negative 
energy, I added to V w (p, z) the hard core 



V Na ci(r) 



,.ii 



47re n e„, 12r 12 ' 



(36) 



If we keep only the 1 Jr term of V w (p, z), then the poten- 
tial V w + VnciCI is proportional to 



,11 



Ylr 1 



(37) 



which has a minimum at r = tq. I took this parameter 
to be ro = 0.51 nm which is the location of both the 
outer maximum of the NaCl-in- water correlation function 
g(r) and also the outer minimum of the (scpism plus 
SIf) potential energy of Na + -Cl~ in water [50 . This 
choice of ro allows the Na + and CI - ions to keep their 
hydration shells; 97% of them do keep their hydration 
shells at 100 mM and 25 C [50; . To prevent the chloride 
ions from falling into the positive charges of the arginines, 
I added a similar term to the R-Cl potential but used 
the somewhat larger value of ro = 0.7 nm to account for 
the more spread-out charge of the bidentate guanidinium 
group. 

Since I treated the water and the lipids as dielectrics, 
I took the potential energy of a sodium or chloride ion of 
charge q — ±|e| to be proportional to the transmembrane 
voltage V ce u reduced by the ratio of the two permittivities 
and by the ratio of the height z to the thickness t of the 
lipid bilayer 



V(z)=q^-jV ceU . 



(38) 



This energy is small compared to kTfj. Even for V ce u = 
100 mV and z = t, it's only fcTb/10. I used the nominal 
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FIG. 5: Snapshot of 94 sodium ions (pluses, red), 106 chloride 
ions (stars, blue), and a 12-mer random coil of oligoarginine 
R 12 (x's, green) after 50,000 sweeps. The coordinates are in 
nm. 



value of 60 mV for V ce u in my simulation of the effect of 
the salt on the transmembrane potential. 

The Monte Carlo code measured the transmembrane 
potential 



TABLE III: The voltage differences AV Na ci (mV) due to the 
156 mM Na + and CP ions near an oligoarginine as an 
a-helix, a random coil, or a /3-strand. The resting transmem- 
brane potential AV C£ u is not included, nor that AVcpp due 
to the oligoarginine. 



TV 


a-helix 


R^ random coil 


R N /3-strand 


5 


168 ±4 


148 ±3 


143 ±4 


6 


202 ±4 


189 ±4 


178 ±1 


7 


233 ±4 


218 ± 3 


204 ± 2 


8 


270 ±3 


244 ± 1 


226 ± 5 


9 


306 ±6 


275 ±2 


250 ±4 


10 


339 ±2 


297 ±5 


266 ±4 


11 


370 ±5 


327 ±3 


281 ±4 


12 


406 ±4 


353 ±4 


303 ±2 



TABLE IV: The voltage differences AVcpp + AV Na ci (mV) 
across the plasma membrane induced by an R^ oligoarginine 
as an a-helix, a random coil, or a /3-strand and by the ions of 
156 mM Na + and Cl~ reacting to it. The resting transmem- 
brane potential AV ce ii is not included. 



N 


R^ a-helix 


R^ random coil 


R^ /3-strand 


5 


- 144 ± 4 


- 154 ±3 


- 148 ± 4 


6 


- 174 ± 4 


- 173 ± 4 


- 168 ± 1 


7 


- 206 ± 4 


- 201 ± 3 


- 189 ± 2 


8 


- 232 ± 3 


- 228 ± 1 


- 199 ± 5 


9 


- 256 ± 6 


- 246 ± 2 


- 205 ± 4 


10 


-281 ±2 


- 260 ± 5 


- 210 ±4 


11 


- 306 ± 5 


- 260 ± 3 


-218 ±4 


12 


- 323 ± 4 


- 261 ± 4 


- 213 ±2 



AV NaC i = V NaC i(0, 0, -t) - V NaC i(0, 0, 0) (39) 



F. Summary of the Model 



due to the salt ions of the active volume. It used the first 



100 terms of the potential V c (p, z) in the cytosol (25) to 
compute Vjv a c/(0, 0, — t), and it used the first 100 terms 



of the potential V w (p, z) in the extracellular medium (24) 
to compute VjsiaCii^i 0,0). The errors of truncation were 
negligible. 

Each run started by assigning random positions to the 
94 Na + ions and the (94 + n) Cl~ ions of the active 
volume and to the 470 sodium and 470 chloride ions of 
the passive volume. After this initialization, the code did 
25,000 thermalizing sweeps in which every Na + and CP 
ion of the active volume was allowed to move as much as 
l/4th of its range in each direction. After thermalization, 
the code measured the transmembrane potential every 10 
sweeps for a total of 2500 measurements. Five runs were 
done for each number 5 < n < 12 of arginines. The 
resulting transmembrane potentials AVnuci due to the 
salt are listed in mV in Table IllTl 



The code took snapshots of the distributions of the 
sodium and chloride ions every 2500 sweeps after ther- 
malization. Fig. [5]displays the last snapshot (after 50,000 
sweeps) of 94 NaT, 106 CP, and 12 Rs in a random coil. 
The coordinates are in nm. 



In the present model of CPP transduction, the trans- 
membrane potential is the sum of three terms — the rest- 
ing transmembrane potential AV ce u of the cell in the 
absence of CPPs, the transmembrane potential AVcpp 
due to an oligoarginine, and the transmembrane poten- 
tial AVjvaC/ due to the counterions of the extracellular 
medium 



Ay = ap cc 



AV C 



CPP 



AVi 



NaCl- 



(40) 



The dominant term is the one AVcpp due to the CPP; it 
is nearly twice as big as the one AVnuCI due to the salt 
and of opposite sign. The resting transmembrane poten- 
tial AVceii of the cell, which arises mostly from the phos- 
phatidylserines of the inner leaflet, augments the sum 
AVcpp + AV~N a ci by some 10-50% depending upon the 
CPP's charge and the value of AV ceU . This salty CPP- 
PS capacitor increases the transmembrane potential V 
and so elevates the Boltzmann factor e~ B ( r "( kT ' and so 
increases the probability of pore formation — at least for 
R^s with enough arginines. It is hard to be quantitative 
here because the voltage required to form a pore depends 
upon the duration of the voltage, the radius of the pore, 
and any defects or fluctuations in the membrane. 
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In its use of an electric field and of the binding of the 
CPPs to the phosphate groups of the phospholipids of 
the outer leaflet, the model has something in common 
with the adaptive-translocation model of Rothbard, Jes- 
sop, and Wender (23]; in its invocation of electroporation, 
it has some overlap with the work of Binder and Lind- 
blom 51J; in its use of neutral dipolar PC & SM head 
groups it is somewhat similar to the work of Herce and 
Garcia [S2] and of Tang, Waring, and Hong ^5 . The 
key distinctive feature of the present model is its under- 
pinning of continuum electrostatics and its quantitative 
synthesis of the contributions of the CPP, the salt, and 
the phosphatidylserines which combine to form a salty 
CPP-PS capacitor with a voltage high enough to cause 
reversible electroporation. 



V. COMPARISON WITH EXPERIMENT 
A. Empirical Upper Limit on Size of Cargo 

Various groups have found that cell-penetrating pep- 
tides cannot transduce cargos of more than about 50 
amino acids |18j . an upper limit that surely varies with 
the cell, the CPP, and the cargo. In the transduction ex- 
periments (i 13 El EllElig[n3ini 1131131111113 IIS] aimed 
at eventual therapies, the heaviest cargo was 33 amino 
acids. The present model based on molecular electro- 
poration offers a qualitative explanation for this upper 
limit. 

The masses of larger cargoes of A = 50-100 amino 
acids of 130 Da each together with N arginines have 
masses Mm, a of from Mn : a — 0.15627V + 6.5 to Mm, a — 
0.15627V + 13 kDa. Our previous formula Q gives the 
lower bounds on the radii of such proteins that run from 
1.29 to 1.58 nm for CPPs of N = 10 arginines. But the 
energy E(r) of a pore rises with its radius r as shown by 
Eq.(|19| and so the chance of pore formation falls with 
the pore radius as shown by Fig. [l] So the chance of a 
pore forming that is big enough for a cargo much larger 
than 50 aa is small. Such cargoes can't easily fit through 
the pores that are most likely to form. 



B. Experiments with Mouse Myoblasts 

Tiinncmann et at [SI used confocal laser scanning mi- 
croscopy to measure the ability of the L- and D-isoforms 
of oligolysine and of oligoarginine to carry fluorophores 
of ~ 400 Da into live C2C12 mouse myoblasts within 
one hour. They found that oligoarginines transduced 
the fluorophores much better than oligolysines and that 
more arginines meant faster transduction, with L-R9 and 
L-R10 doing better than their shorter counterparts as 
shown in Table \V\ They also found that the D-isoforms 
worked better than the L-isoforms and that transduc- 
tion rose with the CPP concentration faster than linearly, 
which may suggest a cooperative effect. 



TABLE V: The fractions of mouse C2C12 myoblasts trans- 
duced by oligoarginines L-R^ at three concentrations (/iM). 



N 


10 fiM 


50 fjM 


100 fjM 


5 


0.0 ±0.01 


0.0 ±0.01 


0.01 ±0.01 


6 


0.0 ±0.01 


0.04 ±0.03 


0.28 ±0.03 


7 


0.0 ±0.01 


0.18 ±0.03 


0.75 ±0.03 


8 


0.02 ±0.02 


0.31 ±0.06 


0.85 ±0.04 


9 


0.05 ±0.04 


0.42 ±0.10 


0.90 ±0.04 


10 


0.70 ± 0.04 


0.73 ±0.04 


0.91 ±0.04 


11 


0.81 ±0.07 


0.83 ±0.03 


0.92 ±0.04 


12 


0.90 ±0.04 


0.92 ±0.02 


0.99 ±0.02 



TABLE VI: The voltage differences AV (mV) across the 
plasma membrane induced by two oligoarginines separated 
by 2 nm and configured as a-helices or as /3-strands. The 
solution was 156 mM NaCl as in Table |TV| Includes a resting 
potential AV ce u ~ —30 mV. 



N 


2 a-helices 


2 /3-strands 


5 


- 249 ± 3 


- 238 ± 4 


6 


- 286 ± 7 


- 271 ± 6 


7 


- 333 ± 5 


- 284 ± 4 


8 


- 365 ± 6 


- 295 ± 3 


9 


- 401 ± 11 


- 304 ± 4 


10 


- 447 ± 3 


- 308 ± 7 


11 


- 483 ± 2 


-310 ±6 


12 


- 509 ± 5 


- 308 ± 5 



The present model is consistent with these experimen- 
tal facts and explains them as follows: The oligoarginines 
crossed cell membranes more easily than the oligolysines 
because they were better able to bind to the phosphate 
groups of the PCs and SMs in the outer leaflet; the 
oligolysines were not able to form a stable upper plate of a 
salty CPP-PS capacitor. CPPs with more arginines were 
transduced more rapidly because with more arginines 
they could bind to more PCs and SMs and because their 
higher charges led to higher transmembrane potentials, 
as noted in Table [TTJ The D-isoforms worked better than 
the L-isoforms because the capacitor mechanism is in- 
sensitive to the chirality of the amino acids and because 
proteases were less able to cut them. To check for a 
cooperative effect, I ran some Monte Carlo simulations 
in which two oligoarginines were as close as 2nm. In 
these simulations, I set ro = 0.55 nm for NaCl and 0.7 
nm for R-Gdm. The resulting transmembrane potentials 
AV are listed in Table |VT| for AV oeU = -30 mV. They 
are higher than those due to a single R™, which appear 
in Table IV (even after AV ce ii = —30 mV is added in). 
Thus higher CPP concentrations accelerate transduction 
because they increase the odds of two or more CPPs 
attaching to nearly the same spot on the outer leaflet. 
There is also the possibility that under physiological con- 
ditions two oligoarginines might form an anti-parallel (3- 
sheet [53] . Such /3-sheets would entail a cooperative ef- 
fect. 
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This consistency of the capacitor model and its simplic- 
ity lends it some plausibility. But evolution finds what 
works, not what fits neatly into a model, and so other 
CPPs with different cargos may enter different cells by 
different mechanisms. In particular, this model may not 
apply to model amphiphilic peptides (MAPs). 

VI. THREE TESTS OF THE MODEL 

One way to test the model would be to compare the 
rates of polyarginine transduction in wild-type cells and 
in those that have little or no phosphatidylserine (PS) 
in their plasma membranes. If PS plays a role as in the 
model of this paper and augments the transmembrane 
potential by 10-50%, then the transduction of polyargi- 
nine fused to a cargo of less than 30 amino acids should be 
somewhat faster in the wild-type cells than in those with- 
out PS in their plasma membranes. Mammalian cell lines 
that are deficient in the synthesis of phosphatidylserine 
do exist [551 [55J [57J EH ES], but they appear to have 
normal levels of PS in their plasma membranes |59j— 
presumably due to a lower rate of PS degradation [BP] , 

Another test would be to construct artificial asym- 
metric bilayers [HTJ [521 [53] with and without PS on the 
"cytosolic" side and to compare the rates of CPP-cargo 
transduction. If the present model is right, then the 
rate of transduction should be somewhat higher through 
membranes with PS on the cytosolic side than through 
membranes with no PS or with PS on both sides. 

If CPPs do enter cells via molecular electroporation, 
then it may be possible to observe the formation of tran- 
sient (< 1 ms) pores by detecting changes in the conduc- 
tance of the membrane [45] . Such measurements would 
be a key test of the model and, if done on an artificial 
membrane, would let one determine both whether CPP- 
transduction is related to the presence of PS on the cy- 
tosolic side of the membrane and whether it proceeds via 
molecular electroporation. 



VII. SUMMARY 

Cell-penetrating peptides (CPPs) can carry into cells 
cargoes with molecular weights of as much as 3,000 Da — 
much greater than the nominal limit 500 of the "rule of 
5" [53]. Therapeutic applications with well-chosen pep- 
tide cargoes of 8-33 amino acids are described in refer- 
ences [l[5J[51[7J[Hll[ini[IIl[Il[I3[Il[I3[I5]. 

Sec. II VI describes a model in which molecular electro- 
poration and phosphatidylserines (PSs) play key roles in 
the transduction of CPP-cargo molecules. In this model, 
one or more positively charged CPPs on the outer leaflet 
and the negatively charged PSs under it on the inner 
leaflet form a kind of capacitor with a transmembrane 
potential in excess of 180 mV for a single CPP of nine 
arginines. This transmembrane potential increases the 



chance of the formation of electropores through which the 
CPP and its cargo can enter the cell. The model is consis- 
tent with the empirical upper limit on the cargo of about 
50 amino acids and with data |31j on how the probabil- 
ity of transduction of polyarginine CPPs into mouse my- 
oblasts depends upon the concentration of the CPP-cargo 
molecules and the number of arginines in each CPP. 

The model predicts that mammalian cells that lack 
phosphatidylserine in their plasma membranes transduce 
polycations less well than those that do, that artificial 
asymmetric bilayers with PS on the cytosolic side trans- 
duce polycations better than ones without PS, and that 
the passage of CPPs should be accompanied by transient 
rises in the conductance of the membrane of the cell or 
BLM. 



APPENDIX A: FIRST ELECTROSTATIC 
PROBLEM 

Here we derive in the continuum limit the electrostatic 
potential V(p, (f>, z) in cylindrical coordinates due to a 
charge q on the z-axis at the point (0,0, h) at a height h 
in the extracellular environment above the phospholipid 
bilayer of a eukaryotic cell for the case in which the height 
h does not exceed the thickness t of the lipid bilayer. 

In electrostatic problems, Maxwell's equations reduce 
to Gauss's law 

VD = i o (Al) 

which relates the divergence of the electric displacement 
D to the density p of free charges (charges that are free 
to move in or out of the dielectric medium — as opposed 
to those that are part of the medium and bound to it by 
molecular forces), and the static form of Faraday's law 

V X E = (A2) 

which implies that the electric field E is the gradient of 
an electrostatic potential 

E = — VV. (A3) 

Across an interface with normal vector n between two 
dielectrics, the tangential component of the electric field 
is continuous 

n X (E 2 - Ei) = (A4) 

while the normal component of the electric displacement 
jumps by the surface density a of free charge 

n • (D a - Di) = a. (A5) 

In a linear dielectric, the electric displacement D is 
proportional to the electric field E 

D = eE (A6) 
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and the coefficient e is the permittivity of the material. 
The permittivity e(m) of a material m differs from that 
of the vacuum e by the electric susceptibility \ and by 
the relative permittivity e m 



e(m) = £o + X = £m e . 



(A7) 



The relative permittivity e m often is denoted K m . 

The lipid bilayer is taken to be flat and of a thickness 
t ~ 5 ran. The relative permittivity of the lipid bilayer is 
eg ~ 2, that of the extra-cellular environment is e w ~ 80, 
and that of the cytosol is e c 80. 

We use the method of image charges. The charge 
q at (0,0, h) will generate image charges at the points 
r — (0, 0, 2nt ± h) in which n runs over all the integers. 
The cylindrical symmetry of the problem ensures that the 
potential is independent of the azimuthal angle 4> and so 
can depend only upon p and z. With p 2 = x 2 + y 2 , the 
potential in the lipid bilayer is 



Vi(p,z) 



1 



9o_ 

vV + (* 



hf 



(A8) 



E E 

s— ±1 O^n— — oo 



9n,. 



^/p 2 + (z- {2nt + sh)f 



while that in the extracellular environment is 



V w (p,z) 



1 







47rene 



E E 



s— ±1 n- 



^ Vp 2 + (^~(2^ + ^)) 2 
(A9) 



and that in the cytosol is 
1 



V c {p,z) 



fc c 



EE 



s— ±1 n- 



;J Vp 2 + (^~ (2ni + s/i)) 2 ' 
(A10) 

The continuity (|A4|) of the transverse electric field £" p 



and that ( A5 ) of the normal displacement D z across the 
planes z — and z = — t imply that the coefficients qo, 
Qn,si Qwn.s, and q C n,s must satisfy for n > and s = ±1 
the relations 



9n,s "T 9— n, — s 


Qk; — n. — s 


(All) 


9n,s 9— ro,— s 


Qw —n,—s 


(A12) 


9«,s + g_( n+ i),_ s = 


Qcn.s 

e c 


(A13) 


9n,s ^ 9-(n+l),-s = 


Qcn , s 


(A14) 



as well as the special cases 



the last two of which imply that 

5-1,1 = 9c0 ,-i = 0. 



(A21) 



The four equations (All - A14l tell us that for n > and 
s = ±1 



yn,s — 


Qw—n,- 




+ H 


Q—n. — s 


~ Qw—n, — s 








e c + eg 


Qn.s — 


n Qcn.s 




£ c ~ Q 

/7„ „ 



from which we can infer that for n > 

yVi,s P9—n,—s 
and that for n > 1 and s = ±1 
J\n— 1 



9n,s = (pp)" 91 

(l+p')^')^ 1 



Qcn , s 
Qw—n,s 



91, 

n— IJn— 1 



-(i+p)p*-V 



9i, 



The four relations (A15 A19 1 imply that 

9u>0,-l = P?MlO,l 

9o = (1 - p)9™o,i 
9c0,i = (1 -P)(l+P')9w0,i 
9-1,-1 = -(1 -P)p'9w0,i- 



(A22) 
(A23) 
(A24) 
(A25) 

(A26) 

(A27) 
(A28) 
(A29) 
(A30) 

(A31) 
(A32) 
(A33) 
(A34) 



Gauss's law (|Al| applied to a tiny sphere about the phys- 
q w o.i = 9- (A35) 



ical charge q gives 



This identification and the four equations (A31 A34) tell 
us that 



9wo,-i = pq 

qo = (i -p)q 

q c0 ,i = (l-p)(l+p')9 

9-1,-1 = -(i -p)p'q- 



(A36) 
(A37) 
(A38) 
(A39) 



Equations ( A26 - A39 ) allow us to relate all the coeffi- 



cients for n > to q w o,i = q and to <7i,-i = 9-i,i = 0: 











9n,l 


9™ o,i 


+ 9m0,-i 


= —qo 

eg 


(A15) 


9n,-l 


9™o,i 


— 9m0,-i 


= 9o 


(A16) 


9-n,l 


9o 
9o 


+ 9-i,-i 
- 9-i,-i 


— — 9e0,l 

e c 

= 9c0,i 


(A17) 
(A18) 


9-n,-l 
Qw — n , 1 




9-1,1 


eg 

— — 9c0,-i 
e c 


(A19) 


Qw —71,-1 

Qcn,l 




9-1,1 


= — 9c0,-i 


(A20) 


Qcn,-1 



( PP T(1-P)q 
(pp')"-^!,-! - 

-p"-y«-V,-i = o 
-p»-y»(i- p ) y 

-(1 +p)p"-V" _1 9i,-i = (A44) 
-(1 -p 2 y i -V™g (A45) 
(l-p)(l+p')(p P Tq (A46) 

(l+j>')(pp') n - 1 gi ) -i = 0. (A47) 



(A40) 
(A41) 
(A42) 
(A43) 
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The electric potential due to a charge q in the extra- 
cellular environment a distance h above a lipid bilayer of 
thickness t then is 



where e wc — (e w + e c )/2. With this approximation and 
with e w w e c ~ 80, the potential (Bl I is 



Vt(p,z) 



q 



i 



47re e„,£ ' \ y/p 2 + {z- 2nt - hf 



n=0 



^p 2 + (z + 2(n+l)t + h) 2 



(A48) 



in the lipid bilayer. That in the extra-cellular environ- 
ment is 



V w (p,z) 



1 M 



P 



ATre a e w \r p 2 + (z + hf 



(A49) 



in which r is the distance from the charge q. Finally, the 
potential in the cytosol is 



V c (p,z) 



Aire 



■E 



(ppT 



o£we£e c ^-~' Q \J p 2 + (z — 2nt — h) 2 



(A50) 



oo p2n 



V a {r k ) » cr / dp / 



OO /*27T 



(B3) 



47re £u, „/r Jo V/° 2 + r fc ~ 2 / 3r fe cos 



(fp / ( — ) -Pn(cOS( 



4^606^, J R Jo ^ \P 

The ro = term in this sum is an infinite constant, which 
we drop because it does not affect the containment of the 
PSs within the disk. The remaining terms are 



V a {r k ) 



Aire, 
aR 



poo oo 

eoe«, Jr ^ 



2n 



2n 



1 2 



\ " 



2e e w ^— ' 2n - 1 

n— 1 



(2n) 



(n!) 2 2 2 ™ 



(I) 2 ". (B4, 



APPENDIX B: SECOND ELECTROSTATIC 
PROBLEM 



Here I approximate the electrostatic potential V a (rk) 
within a disk of radius R due to a uniform charge density 
a of phosphatidylserines (PSs) outside the disk. 

The negative charges of the PSs are taken to lie on 
the interface between the cytosol and the lipid bilayer. 
The role of this potential V a is only to keep the mutual 
repulsion of the PSs inside the disk from driving them 
too much toward the perimeter of the disk. So an exact 
expression for V a is not needed. Any formula for it will 
involve an integral of a over distances that run to infinity. 
My approximation is to set the thickness t of the bilayer 
equal to zero. In this limit, the effective potential felt by 
a PS at r k is 



/>00 />27T 

V a (r k ) = a / dp / d<t>"V ws {\p-r k \) 
Jr Jo Q 



(Bl) 



in which p = p(cos <fi, sin 6, 0), and V ws (\p — r k \) is the 
potential V w (p, z) of Eq. (|A49[) for z = h = t = 



V ws {\p - r k \) 



1 



47re e wc |p- r k \ 



(B2) 
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